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Single photon state generation from a continuous wave non-degenerate optical 
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We present a theoretical treatment of conditional preparation of one-photon states from a con- 
tinuous wave non-degenerate optical parametric oscillator. We obtain an analytical expression for 
the output state Wigner function, and we maximize the one-photon state fidelity by varying the 
temporal mode function of the output state. We show that a higher production rate of high fidelity 
Fock states is obtained if we condition the outcome on dark intervals around trigger photo detection 
events. 
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I. INTRODUCTION 

Light has come to play an important role in quantum 
information sciences, and this puts focus on the ability 
to prepare and manipulate state vectors of light, i.e., to 
restrict the field to a few modes or a single field mode, 
and to interact exclusively with the selected modes. In 
optical and microwave cavities, many experiments have 
already been done on pure quantum state manipulation 
in isolated field modes interacting with single atoms, but 
there is also an obvious interest to produce traveling light 
fields with controllable pure state properties. 

It is not easy to manipulate quantum states of light, 
and one very successful strategy has thus been to pro- 
duce a certain "easy" state, and to use a quantum mea- 
surement to project this state on the desired quantum 
state. Following the theoretical proposals by Dakna et 
al [l[ , see also 0, H| , conditioned on appropriate detec- 
tion events, single- and two-mode squeezed states of light 
can be transformed into Fock states and Schrodinger cat- 
like states, as demonstrated with good fidelity with single 
mode pulses of light in [J, 0, @, Q- I.e., a large number 
of pulses are generated, and measurements, carried out 
on each of them with some success probability, ascertain 
the preparation of the desired state. 

A scheme to produce single photon states and other 
single mode quantum states from a continuous wave field 
was demonstrated in experiments [8|, where a pho- 
ton counter registered the intensity of a small fraction of 
a beam of squeezed light, causing the remaining beam 
to have a high single photon amplitude in a well local- 
ized mode. A theoretical analysis of that experiment 
was given in llCfl. In the present paper we generalize 
the approach of [10| to protocols involving twin beams of 
light generated from a non-degenerate optical paramet- 
ric oscillator (OPO). The non-degenerate OPO converts 
photons from a pump beam into pairs of photons in a 
pair of output modes, and a conditional detection of a 
single photon in one beam (which we shall denote the 
trigger beam) results in the presence of a single photon 
in the other beam (the signal beam). A possible exper- 
imental setup is sketched in figure [TJ In contrast to the 



pulsed case where the field mode occupied is governed by 
the pulse shape, in the continuous wave case, the precise 
temporal modes occupied by the one-photon states in the 
signal beam have to be specified. The largest one-photon 
state fidelities are obtained by choosing the modes opti- 
mally, and we present results of a variational procedure 
and of a numerical optimization. 

In Sec. II, we introduce the quantum correlation func- 
tions of the twin beams of interest, and we show how 
the (Schrodinger picture) quantum state populating an 
arbitrary field mode can be readily obtained from the 
(Heisenbcrg picture) field correlation functions. In Sec. 

III, we present analytical results for the phase space 
Wigner function for the signal state conditioned on a 
single photo detection event and for the fidelity and rate 
with which the one-photon states are produced. In Sec. 

IV, we utilize a variational method to optimize the signal 
mode function. Finally, in Sec. V, we analyze the signal 
state conditioned on a photo detection event surrounded 
by an interval with no photo detection events. Sec. VI 
concludes the paper. 



II. TWO-TIME QUANTUM CORRELATION 
FUNCTIONS AND PHASE SPACE WIGNER 
FUNCTIONS FOR GAUSSIAN STATES 

The OPO Hamiltonian is quadratic in field annihila- 
tion and creation operators, and the solution of the time 
evolution is most easily accomplished by a linear (Bo- 
goliubov) transformation of each pair of incident field 
operators at all frequency components. The state there- 
fore retains the Gaussian state character of the vacuum 
state, i.e., the joint probability distributions of all the 
quadrature variables is Gaussian. The state separates in 
pairs of frequency modes for which the quantum state 
can be written down explicitly. Since we shall deal with 
photo detection experiments, we need a representation of 
the field in time domain, and although it is still a Gaus- 
sian state, it is now much more complicated, as the field 
operators at different times are non-trivially correlated. 

In so-called Type I and Type II OPO's, photons are 
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FIG. 1: Experimental setup for one-photon state production 
using a type II OPO. A polarizing beam splitter (PBS) sep- 
arates the source output into two twin beams with opposite 
polarizations. One beam is used as a trigger, and when the 
trigger detector (APD) registers a photon, the generated state 
is investigated by homodyne detection with a local oscillator 
(LO). 



created pairwise in two different frequency components 
and in two different polarization components, respec- 
tively. It is also possible to create photons pairwise in 
beams propagating in different directions. Alternatively 
one can combine the single beam squeezed light outputs 
from two degenerate OPO's on a beam splitter, which 
also produce twin correlated beams, if the OPO outputs 
are squeezed 7r/2 out of phase. All of these approaches, 
if restricted to correlations in only one degree of freedom 
(polarization, frequency, direction) , lead to quantum cor- 
relations described by the same formal correlation func- 
tions for the annihilation and creation operators for the 
two components, denoted in the following by + and — , 
a+(t), a\(t), Cb-(t), and a_(t) , an d are in quite precise 
models of the OPO's given by [lil ] 
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The parameters in these expressions are the non-linear 
gain coefficient e of the OPO, and the decay rate 7 of light 
in the OPO cavity due to leakage through the output 
mirror. It is assumed that there is no loss through the 
other mirror. For fixed 7 we note that e/7 is closely 
related to the mean twin beam intensity (6^.(t)o±(t)) = 
2( £ / 7 ) 2 7 /(l-4(6/7) 2 )). 

The modes relevant to the experiment are the mode 
in which the trigger detection takes place, and 



the mode fait) occupied by the produced state, while 
all other modes are unobserved. We assume the trigger 
detection to take place on a timescale much shorter than 
7 _1 . The precise shape of the detected temporal mode 
is then irrelevant, and we assume 
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where At c is a short time interval and t c is the click time 
for the trigger detection. As discussed in the introduc- 
tion, the mode function f2(t) can be chosen arbitrarily 
(under the constraint J \f2(t)\ 2 dt = 1), and in particu- 
lar it can be optimized to achieve maximal fidelity. The 
optimal mode function is centered around the trigger de- 
tector click time, and from fTJ it is to be expected that 
the temporal extent is of order 7 -1 , which is the time 
that a photon can spend in the cavity and hence be sep- 
arated from its partner in the output fields. We will re- 
turn to this point in sections III and IV, but by now we 
leave f 2 {t) unspecified. The problem thus reduces to the 
one of characterizing the correlations of the single mode 
operators 
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where r] t is the trigger detector efficiency, rj s is the signal 
detector efficiency, and a± lVac are field operators acting 
on vacuum, included to ensure appropriate commutator 
relations of &i and dj. 

Since the state is Gaussian, we immediately know the 
multi-component Wigner function of the trigger and sig- 
nal mode before conditioning on the trigger detector click 
event. Defining the column vector 

{x 1 ,pi,X 2 ,P2) T , 
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of quadrature variables 
Wigner function is 

W v (y) ' 
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(cii = (£i + ipi)/V2) ■, this 



= exp (-y T V 1 y) . 



(7) 



7rVdet(F) 

Wv(y) is parameterized by the covariance matrix V with 
elements Vij = {yijjj) + {ViVi), which are computed ex- 
plicitly by use of equations (HI) and ((SJ), the time de- 
pendent mode functions fi(t), and the two-time cor- 
relation functions (J7XJ) - Since it follows from (JT|) that 
(diSi) = {0,20,2) = {a\a,2) = 0, we have 

2(a{a 1 ), (8) 

(9) 

= V21 = V 34 = V 43 = 0, (10) 
= V 31 = -V24 = -V i2 = 2Rc((a 1 a 2 )), (11) 
V u = V 23 = V 32 = 2Im«aio 2 )). ( 12 ) 
We note that we use the same symbols for quadrature op- 
erators and for the corresponding real variable arguments 
in the Wigner function. 
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III. ANALYSIS OF THE STATE CONDITIONED 
ON A SINGLE CLICK EVENT 

We first determine the output state Wigner function 
conditioned on a single trigger detector photo detection 
event. We shall assume that the trigger detector pro- 
duces click events with a rate proportional to the mean 
photon number, but it does not by its outcome distin- 
guish between different photon number states. A click 
detection of such a detector is accompanied by the ap- 
plication of the corresponding field annihilation operator 
on the quantum state of the system, i.e., by application 
of the annihilation operator from the left and its adjoint 
creation operator from the right on the density matrix, 
followed by a partial trace over the otherwise unobserved 
modes. These operations are mapped to certain differ- 
ential operators on the Wigner function followed by 
integration over the corresponding phase space variables. 
This results in the following Wigner function for the con- 
ditioned signal state 
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FIG. 2: Wigner function conditioned on a single trigger de- 
tector click event for e/7 = 0.02, r/ a = 1, and f2(t) = 
\Zt72 exp(— j\t — t c \/2). The one-photon state fidelity is 
0.9921, and the negativity at the origin is W c u c k(0, 0) = 
-0.3133. 
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where N c u c k is a normalization constant. By partial in- 
tegration (|13[) reduces to 



Wdick (x2,Pi) = Ndick J diridpi- (x\ +pl - l) W v {y) 

= (A 1 +A 2 {x 2 2 +pl))e- A ^ +p i\ (14) 

where the coefficients A\, A2, and A3 are given in terms 
of the covariance matrix elements (JHHUJ 
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W c ii c k depends on the signal mode function f2(t), the 
signal detector efficiency r) s , and the non-linear gain co- 
efficient e of the OPO through V, but W c iick does not 
depend on the trigger detector efficiency r\ t . Since we 
condition on a click event, the only effect of an ineffi- 
cient trigger detector is to reduce the production rate. 
Figure [2] shows W c u c k f° r e /l — 0.02, -q s = 1, and 
fz(t) — y/ r y/2exp(— r y\t — t c \/2) (in the next section we 
prove that this is the optimal mode function for e/7 — ► 0). 
The fidelity, by which the conditioned Wigner function 
W c iick resembles the Wigner function for a one-photon 
state 
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and for the Wigner function in figure [5] we find 
= 0-9921. 

For completeness we present in appendix A analytical 
results for the Wigner function and the one-photon state 
fidelity for the degenerate OPO setup analyzed in [Icj |. 

The optimal signal mode function (denoted f p(t)) ob- 
tained from numerical optimization of the one-photon 
state fidelity (fT6|) is shown in figure [3] for i] s — 1 and 
different values of e/7. (The optimization is over all real 
functions, which will be justified in the next section.) 
The tendency is that the width of the mode function de- 
creases, when the beam intensity is increased, and a dip 
appears on each side of the peak. However, the difference 
between the optimal mode function for e/7 = 0.1 and for 
e/7 = is quite small. 

The optimized fidelity is shown in figure [5] for rj s = 1 
and ?7 S = 0.8 (solid lines). The fidelity decreases with 
increasing e because the two-photon state contribution to 
the output state increases when the intensity increases. 
The effect of reducing rj s from 1 to 0.8 is roughly to reduce 
the fidelity by 20 % in the small e/7 region. For larger 
e/7 the reduction is smaller because the lower efficiency 
reduces the intensity seen by the detector. The dashed 
curves in the figure show the fidelity calculated for the 
mode function f 2 (t) = y / "f/2exp(— j\t — t c \/2), and it 
is apparent that this analytical approximation is close 
to optimal for e/7 < 0.1. This region is also the most 
interesting region since the fidelity is low for larger e/7. 

The production rate r (i.e. the trigger detector click 
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FIG. 3: Optimal mode function f op (t) for e/7 = (solid 
line), e/7 = 0.1 (dashed line), e/7 = 0.2 (dot-dashed line), 
and e/7 = 0.3 (dotted line). The signal detector efficiency is 
r]s = 1. 




FIG. 4: Fidelity calculated from the numerically optimized 
mode function (solid lines) and from the e = mode function 
in equation (|20|l (dashed lines). The curves approaching 1 for 
small e are for perfect detection rj s — 1, and r) s = 0.8 for the 
curves approaching 0.8 for small e 



rate) is r = P c iick/ At c , where P c u c k is the probability to 
observe a click in the trigger mode function. We assume 
that (ajai) << 1, which is valid if At c is small compared 
to the mean temporal distance between the photons in 
the trigger beam. The probability to observe a click in 
the trigger mode is then the expectation value of the 
number of photons in that mode, and hence 
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Figure [S] shows the production rate as a function of e/7 
for 7]t = 1. The rate decreases when e/7 decreases and 
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FIG. 5: Production rate r versus e/7 for r\t = 1. 



is zero when e/7 = 0, and hence, to achieve an accept- 
able production rate, e/7 needs to be sufficiently large. 
However, according to figure [4] this limits the fidelity. In 
Sec. V we analyze the possibility to reduce this problem 
by conditioning on clicks surrounded by intervals with 
no click detection events, but first we consider the mode 
function optimization in greater detail. 



IV. VARIATIONAL OPTIMIZATION OF THE 
SIGNAL MODE FUNCTION 

In this section we optimize the signal mode function by 
analytical variation of the one-photon state fidelity (|16[) . 
We write the signal mode function as /2(f) = \f2(t)\e l6 ^\ 
where 9{t) is real. |/2(*)| and 6(t) may be varied inde- 
pendently, and since the constraint J |/ 2 (i)| 2 dt = 1 does 
not involve 9(t), the optimal choice of 6(t) fulfils that the 
variation of F\ vanishes, when 9(t) is varied, which leads 
to 

"\f2(t)\\h(t')\ S in(6(t>)-6(t)) ( Cl (f 2 ){aUt)a-(t')) 
+ c 2 {f 2 )(a+(t c )a. (t)) (S+ (t c )fi_ (*')>) df = 0, (18) 
where 
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Note that c\ and C2 depend on the signal mode function 
through the covariance matrix elements, but they do not 
depend on time. Note also that V\\ — 1 and V 2 3 + V 2 4 are 
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both proportional to r) t At c , while V33 is independent of 
rj t and At c . 

For a general \f 2 (t)\ the solution to (fT8|) is 6{t) = 
constant + Ntt, where N is an integer that may be dif- 
ferent at different times. (| 18[) is also fulfilled in the limit 
where 6(t) varies infinitely fast with time, but this leads 
to V33 = 1 and Vi 3 + V 2 A = 0, and hence W e u c k reduces 
to the Wigner function for vacuum and F\ = 0. Thus the 
phase is constant except for possible sign shifts. Adding 
a time independent constant to 9(t) is equivalent to ro- 
tating the signal state Wigner function in phase space, 
but since the desired \n — 1) signal state Wigner func- 
tion depends on x\ and p\ through x\ + p 2 , it is rota- 
tionally symmetric, and hence the fidelity is unchanged 
by this transformation. Consequently, we can choose 
fop(t) = ±\f op (t)\, and it is thus sufficient to optimize 
over all real functions. 

Restricting fi{t) to be real, the constraint is 
J f2(t) 2 dt — 1, and this is taken into account by introduc- 
ing a Lagrange multiplier £ and demanding the variation 
of Fi — £ j f2{t) 2 dt to vanish when f2(t) is varied. This 
leads to the following integral equation for the optimal 
signal mode function 
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For e = the first term in (TT5]) vanishes, and we have 



lim fop(t) = exp (~\t - t c 



(20) 



(|20[) is the optimal zero intensity signal mode function, 
and for 77,, = 1 the fidelity is lim e _>o Fi(f op (t)) — 1. For 
e ^ 0, (|2"D|) can be used as the starting point for a nu- 
merical iteration of (|19p . To increase the range of e/7 
values, where the iteration procedure converges, we use 
the mode function N a ((l — a)(f2(t))i+a(f2(t))i-i) in the 
(i + l) th iteration step, where (/2(i))i is the mode func- 
tion obtained from the i th step, a is a number between 
and 1, and N a is a normalization constant. Using this 
method we regain the mode functions in figure [3] 



the trigger detector is turned off outside the conditioning 
click time window, or if we want to calculate the aver- 
age fidelity for a large number of produced states where 
no selection of particular click events (for instance those 
with no other clicks nearby) has taken place. However, 
an improvement in fidelity for fixed e/7 can be obtained 
if we only accept clicks that are surrounded by a certain 
time period T with no clicks since this restricts the two- 
photon and higher photon number state contributions to 
the produced state. The fidelity increase depends on the 
trigger detector efficiency. If r) t — 1, it is certain that no 
photons hit the trigger detector during the period T if 
no clicks are registered, and for large T7 the fidelity ap- 
proaches unity irrespective of the value of e/7 (provided 
Ai c 7 is much smaller than the mean temporal distance 
between photons in the trigger beam). In the opposite 
limit, rj t = 0, a vacuum detection (i.e. no click) gives us 
no information about the produced state, and we regain 
the results from Sec. III. 

For fixed e/7 the conditioning on a no click interval 
results in a smaller production rate because some clicks 
are now disregarded, but since the fidelity is increased, 
we can increase e/7 a little while still obtaining a larger 
fidelity than in Sec. Ill, and if T7 is not too large, this 
leads to a larger production rate. In the region where the 
rate is limited by the number of photons produced by the 
source it thus turns out that an increase in production 
rate for fixed fidelity or an increase in fidelity for fixed 
production rate can be obtained. 

In the following we first calculate the one-photon state 
fidelity for the signal state conditioned on one click detec- 
tion and an arbitrary number of no click detections, and 
after that we derive an equation for the corresponding 
production rate. Finally we present and discuss numeri- 
cal results for the fidelity and the production rate. 



A. Fidelity 

To determine the one-photon state fidelity for the pro- 
duced state conditioned on a click event surrounded by 
a time interval with no click events, we replace the con- 
tinuous time argument by a discrete set of box shaped 
temporal trigger mode functions and assume the trigger 
detections to take place in these modes. The m = T/At c 
vacuum detection modes are labeled by the numbers 
3, 4, . . . , m + 2 and are included in the covariance ma- 
trix. The unconditioned Wigner function is then 



V. CONTINUOUS DETECTIONS WITH 
(EFFICIENT) TRIGGER DETECTOR 



W v (y) 



1 
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In Sec. Ill we calculated the one-photon state fidelity 
for the output state conditioned on a single click event. 
Conditioning solely on one click event means that we 
average over all possible outcomes of the trigger detec- 
tor measurements outside the small time window where 
the conditioning click occurs. This method is correct if 



where y = (x\,pi, x 2 ,p 2l ■ ■ ■ , x m+2 ,p m +2) T ■ We first de- 
termine the Wigner function for the state conditioned 
on all the vacuum detections but not the click. A vac- 
uum detection results in a projection of the relevant 
mode on the vacuum state, which in terms of Wigner 
functions corresponds to multiplication with the vacuum 
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state Wigner function W n =o(x,p) = fl" -1 exp(— x 2 — p 2 ) 
followed by integration over the relevant quadrature vari- 
ables and renormalization. Thus the Wigner function 
W vaccon for the state conditioned on the m vacuum de- 
tections is 



accon [p^l i Pi i %2 ? P2 
m+2 



N, 



vaccon 



II / dx t d Pt W n=0 (x llPl )jW v (y), (22) 

i=3 ' 



where N vaccon is a normalization constant. Since the vac- 
uum projection is a gaussian operation, W vaccon is also a 
gaussian function. This means that we can proceed as in 
section III if we just replace the covariance matrix for the 
output mode and the click mode with the covariance ma- 
trix Vvaccon corresponding to the gaussian Wigner func- 
tion W vaccon . If we write the (2m + 4) x (2m + 4) covari- 
ance matrix Va m +4 for the signal mode, the click mode, 
and the vacuum modes as 



2m+A 



v 4 c 



(23) 



where V4 is the 4x4 covariance matrix for the signal mode 
and the click mode, and Vi m is the covariance matrix for 
the vacuum modes, it has been proven [l3| that 



V v 



= v i -c(v 2m + h m )- 1 c T 1 



(24) 



where I^m is the (2m) x (2m) identity matrix. Hence by 
use of 1|24[) we can immediately calculate the improved 
fidelity from (JTHJ) (note that the covariance matrix ele- 
ments that are equal in Sec. Ill are changed by equal 
amounts and are hence still equal). 



B. Production rate 



The production rate is 



P, 
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(25) 



where P vac ,click is the probability to detect no click in 
the m vacuum modes and a click in the click mode. To 
determine P, 



vac.ctic 



lick 



we calculate the overlap of the un- 
conditioned source output state with the vacuum state 
in the m vacuum modes and then the expectation value 
of a\di = \(x\ +Pi — 1) in the click mode 
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xi+pi-i) 



11 W n =o(xi,Pi) W v (y)dy. (26) 



Since we consider trigger modes only the signal mode is 
not included in the (2m + 2) x (2m + 2) covariance matrix 
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FIG. 6: Optimized one-photon fidelity as a function of vac- 
uum detection window T. The trigger detector efficiency is 
rjt — 1 for the dotted line, r\t = 0.4 for the solid line, and 
r/t — for the dashed line, while e/7 = 0.2 for all curves. 



Performing the integration we find 



P 



vac. click 



v /det(/ 2 m+2 + VJ) 

f det ((V^ 1 + J x )red) 

{ det (V x - 1 + J x ) 



- 1 



(27) 



where we have introduced the following notation: Ii m +i 
is the (2m+2) x (2m+2) identity matrix, J is a (2m+2) x 
(2m + 2) matrix with elements J a = 1 for i > 2 and Jij — 
otherwise, V x is the (m+ 1) x (m + 1) matrix consisting 
of all the odd-odd matrix elements of V, J x is an (m + 
1) x (m + 1) matrix with elements (J x )u = 1 for i > 1 



and (J x )i 



otherwise, and red means that the first 



row and the first column of the matrix (corresponding to 
the click mode) have been removed. 



C. Numerical results 

In all the numerical calculations in this subsection rj s = 
1 and At c "/ = 0.02. For e/7 = 0.2 (and rj t = 1) we find 

(0J01) 



2-10 a , and hence the assumption (ajai) << 1 
is valid for the region, where the fidelity is large. Even 
for e/7 = 0.45, (a|ai) = 4- 10~ 2 is still somewhat smaller 
than 1. 

Figure [6] shows the optimized one-photon state fidelity 
as a function of the time interval T in which we de- 
mand no click events to occur. T is chosen symmetrically 
around the trigger detector click time, i.e. we condition 
on no clicks from t = t c - At c /2 - T/2 to t = t c - At c /2 
and from t = t c + Ai c /2 to t = t c + At c /2 + T/2, since 
this gives rise to the largest fidelity increase for a given 
T7. The rate and fidelity expressions are, however, also 
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FIG. 7: Optimal mode function for rjt = 1 (dotted line), ryt = 
0.4 (solid line), and rjt = (dashed line), e/7 = 0.2 and 
T7 = 10. The dot-dashed line is the e/7 = optimal mode 
function given by (|20[) . 
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FIG. 8: Maximal one-photon fidelity versus e for T7 = 10. 
The dotted line corresponds to rjt = 1, while rjt — 0.4 for 
the solid line. The dashed line is for rj t = 0, and it is hence 
identical to the upper (solid) curve in figure [4] 



valid for asymmetrical time intervals (and even for time 
intervals interrupted with periods where the trigger de- 
tector is turned off). As expected the fidelity increases 
when T7 increases if rj t ^ 0, and for rj t = 1 the fidelity 
approaches one for large Tj. For a moderate trigger de- 
tector efficiency rjt = 0.4, the fidelity increase is smaller, 
but still significant. The fidelity increase levels off when 
T7 increases beyond approximately 10 irrespective of the 
value of rjt . The reason for this is that the temporal ex- 
tent of the signal mode function is approximately 10 7 _1 
as is apparent from figure[7l which shows the optimal sig- 
nal mode function for T7 = 10, e/7 = 0.2, and different 
trigger detector efficiencies. For r\ t = the mode func- 
tion is identical to the e/7 = 0.2 mode function in figure 
[31 but when r] t is increased to non-zero values, the mode 
function approaches the optimal e/7 = mode function 
([2TJ]) (the dot-dashed line in the figure). 

The optimized one-photon fidelity as a function of e 
for T7 =10 and different values of trigger detector effi- 
ciency is plotted in figure[51 The dashed curve (rjt — 0) is 
identical to the T — curve in Sec. Ill, so the difference 
between the solid line and the dashed curve shows the 
fidelity increase AFi when T7 increases from to 10 for 
rjt = 0.4. AFi increases with e because the mean tempo- 
ral distance between the photons in the trigger and the 
signal beam decreases when e/7 increases, and hence it 
is more likely to have close clicks. The dotted curve in 
the figure is for perfect trigger detection rjt = 1, and as 
expected it is close to unity for all e/7. 

The production rates corresponding to the solid line 
and the dotted line in figure [5] are shown in figure O and 
the dashed line in figure [9] is the rate for T = and 
rjt = 0.4. The production rate vanishes for zero intensity 
because there are no photons present in the trigger beam 
in that limit. Also in the large intensity limit the rate 
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FIG. 9: Production rate r as a function of e. The dotted line 
corresponds to rjt — 1 and T7 = 10, while the solid line is for 
rjt = 0.4 and T7 = 10. The dashed line is for rjt = 0.4 and 
T7 = 0. 



approaches zero for T^0. The reason is that the large 
photon flux makes it extremely unlikely to observe no 
clicks during some non-zero time interval. It is desirable 
to be on the left side of the maximum of the relevant 
curve in figure [5] because the fidelity is a decreasing func- 
tion of e. Hence for T7 =10 and rj t = 1, e/7 should be 
below 0.20, while for rj t = 0.4, e/7 should be below 0.28. 
The % = 0.4 curve is shifted to the right compared to 
the rjt = 1 curve because a lower detector efficiency corre- 
sponds to a lower photon flux in the beam. The maximal 
rate (r/7 ~ 0.029 for both curves) corresponds to a mean 
temporal distance between accepted clicks of 34 7 _1 . For 
rjt = 0.4 the high fidelity region is e/7 < 0.15. For this 
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region it is apparent from figure [9] that the decrease in 
production rate for fixed e/7 is quite small when T7 is in- 
creased from to 10, while figure [8] shows that the fidelity 
increase AFi is significant. 



VI. CONCLUSION 

In conclusion we have considered conditional single 
photon state generation from a continuous wave non- 
degenerate OPO. We have presented explicit analytical 
expressions for the output state Wigner function, the 
one-photon state fidelity, the production rate, and the 
optimal zero intensity limit temporal mode function for 
the signal state. We have optimized the mode function 
numerically to achieve maximal fidelity, and from a vari- 
ational calculation we determined an integral equation 
for the optimal mode function. By conditioning on clicks 
surrounded by an interval with no clicks, the fidelity was 
increased for fixed e/7, while the production rate was 
decreased slightly. However, accepting a smaller fidelity 
increase by increasing e/7 led to a larger production rate 
if T7 and e/7 were not too large. 

Using our approach, it is possible analytically to cal- 
culate the output state Wigner function conditioned on 
arbitrary sequences of click detections, vacuum (i.e. no 
click) detections, and no detections (i.e. detector turned 
off) in temporal trigger modes of short duration, and the 
corresponding production rate can also be determined 
analytically since both calculations involve integration 
of products of polynomials and Gaussian functions only. 
From the conditioned output state Wigner function it 
is straightforward to compute the fidelity for arbitrary 
states in arbitrary output modes. We are currently work- 
ing on a generalization of our calculations for the produc- 
tion of single photon states to generation of n-photon and 
Schrodinger kitten states. 

The authors acknowledge discussions with Jonas S. 
Neergaard- Nielsen and Eugene S. Polzik, and financial 
support from the European Union through the Integrated 
Project "SCALA". 



APPENDIX A: WIGNER FUNCTION AND 

FIDELITY FOR ONE-PHOTON STATES 
PRODUCED FROM A DEGENERATE OPO 

The setup analyzed in [ic| is obtained from figure [1] 
by replacing the non-degenerate OPO with a degener- 
ate OPO and the polarizing beam splitter with a normal 
beam splitter with low transmission. The two-time cor- 



relation functions for this setup are [Io| 

A 2 - [P- /e^l* - *'! e — A|t— 1'| 




4 \ 2li 
and ([4]) and §5§ are modified to 



+ y/l-Vtbvac(t'))dt', (A3) 



a 2 = j f 2 (t') [y^Ra{t') + - R)a>vac(t') 

+ y/l-VsCvac(t'))dl/, (A4) 

where R is the beam splitter reflectivity, and the opera- 
tors labeled by vac are field operators acting on vacuum. 
As in Sec. IV we can restrict our analysis to real sig- 
nal mode functions. In terms of the covariance matrix 
elements computed from ([AT]) . (TA"2"j) . (|A"3j) . and (TA"4")) the 
Wigner function for the produced state conditioned on a 
single trigger detector click event is 



1 



(C 2 + C 3 xj + dpDe 



-Csx 2 — Cep 2 



(A5) 



where 






C x = 


^(U 33 U44) 5/2 (Ull + 


■V22-2), 


c 2 = 




+ V22 - 2) 


c 3 = 


2^44. 




d = 


2V24V33, 




c 5 = 






c 6 = 







and the one-photon state fidelity is 
2(U 33 U 44 - 1) 



(1 + U 33 )3/ 2 (1 + U 44 ) 3 / 2 
2(2(1 + U 44 ) + 1 - U 33 U 44 )Vi 2 3 

(Vu + V22 - 2)(1 + U 33 )5/2(l + y 44 )3/2 



2(2(1 + U 33 ) + 1 - U 33 U 44 )K 



24 



(Vn + V22 - 2)(1 + U 33 )3/ 2 (l + U 44 ) 5 / 2 



(A6) 
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